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We investigate the phase difference of the sunspot cycles in the two hemispheres and compare it with the latitudinal sunspot 
distribution. If the north-south phase difference exhibits a long-term tendency, it should not be regarded as a stochastic 
phenomenon. We use datasets of historical sunspot records and drawings made by Staudacher, Hamilton, Gimingham, 
Carrington, Sporer, and Greenwich observers, as well as the sunspot activity during the Maunder minimum reconstructed 
by Ribes and Nesme-Ribes. We employ cross-recurrence plots to analyse north-south phase differences. We show that 
during the last 300 years, the persistence of phase-leading in one of the hemispheres exhibits a secular variation. Changes 
from one hemisphere to the other leading in phase were registered near 1928 and 1968 as well as two historical ones 
near 1783 and 1875. A long-term anticorrelation between the hemispheric phase differences in the sunspot cycles and the 
latitudinal distribution of sunspots was traced since 1750. 
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1 Introduction 

In order to deal with the asymmetry of hemispheric sunspot 
distributions, separate observations for the northern and 
southern hemispherical activities are required. Suitable 
monitoring was started by the Royal Greenwich Observa- 
tory in 1874. Now daily information about size, position and 
area of sunspots are comprised by the RGO/USAF/NOAA 
dataset. Unfortunately, both the International Sunspot Num- 
bers and the Group Sunspot Numbers are defined for the 
entire visible solar disk only. 

Applying the cross-recurrence plot technique to the 
smoothed sunspot area data, we compared the time-varying 
phase differences between the northern and southern hemi- 
spheres with the latitudinal distribution of the sunspots 
(Zolotova et al. |2009l l. If the phase difference between the 
northern and southern hemispheres is positive, we refer to 
the northern hemisphere as 'leading in phase'. If the dif- 
ference is negative, the southern hemisphere is 'leading in 
phase' . We observed a long-term persistence of one hemi- 
sphere exhibiting a leading phase, which lasts almost four 
solar cycles, while after that, the opposite hemisphere leads 
in phase. The full period corresponds to the secular cycle 
first described by Gleissberg ('1967V We also verified that 
these long-term variations in the hemispheric phase differ- 
ences are in anticorrelation with the long-term variations of 
the latitudinal distribution of the sunspots - the magnetic 
equator (PuUckinen et al. , 1999V 

Since the analysed time series starts in May 1874, we 
were unable to detect the change in sign of the phase dif- 
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ference before the 12th cycle. In the present paper we ex- 
tend our previous studies of historical hemispheric sunspot 
asynchrony as far back as to the Maunder minimum. Direct 
observations of the Sun by Staudacher (Arlt 2008i |2009al l. 
Hamilton & Gimingham (Arlt I2009bl) . Carrington (fT863T l. 
and Sporer (fT874l [18781 [1880' T886' 1897) are involved. 
Further, we reconstruct and discuss the temporal variation 
of the phase difference and the latitudinal distribution of the 
sunspots during the last 300 years. 



2 Data and methods 

For the analysis we use the following list of daily data: 

- The drawing of the butterfly diagram during the 
Maunder minimum constructed by Ribes & Nesme- 
Ribes ([T993T l. 

- The sunspot drawings by J.C. Staudacher for the pe- 
riod 1749-1796. Their digitization, determination of 
heliographic latitudes and longitudes in the Carring- 
ton rotation frame, and reconstruction of the butter- 
fly diagram was done by Arlt (l2008l |2009al l. These 
data were extended with a few solar observations by 
J. A. Hamilton and W. Gimingham for the period 1795- 
1797 (Arlti2009b). 

- The records by R.C. Carrington for the period 1853- 
1861. His famous book contains the heliographic coor- 
dinates of individual sunspots (Carrington I1863I) . The 
last part of the book represents sunspot group drawings 
for each rotation, starting from the first in 9 November 
1853. One of Carrington's books was presented by the 
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Royal Astronomical Society to the Russian Imperial Li- 
brary (now the Library of the Russian Academy of Sci- 
ence), and one to the Helsinki University Ubrary. 

- The records of heliographic sunspot positions by 
G. Sporerforthe period 1861-1894 (Sporeri87l|T828l 
[18801 [T8861 [18941 . All books consist of a description, 
tables with positions, and drawings of sunspot groups at 
their latitudes; the longitude is not shown in the draw- 
ings, instead, the horizontal axis is used for the date. 
These observations are not ordered by calendar dates, 
but by a sequence of group numbers. This data and 
Carrington's were digitized in the Helsinki Observatory 
about 20-30 years ago and arranged as a function of 
time. 

- The RGO/USAF/NOAA dataset from 1874 to the 
present, originally in the Royal Observatory, Green- 
wich [187097^ . 

Finally, we compiled a list of heliographic positions of 
sunspots for those days, when observations occur (except 
for the Maunder minimum, when only the butterfly dia- 
gram is available; Ribes & Nesme-Ribes I1993I I. To count 
sunspots for each hemisphere we replaced the latitude val- 
ues by plus (+1) and minus (—1) respectively for the north- 
ern and southern hemispheres. Next, we added numbers be- 
longing to the same day but separately for the hemispheres. 

To measure the phases differences of sunspot cycles in 
the two hemispheres we have used the Cross-Recurrence 
Plot (CRP) method (Marwan et al. I2007I I. where two time 
series are embedded in the same phase space. The CRP 
analyses the parallel occurrence of states as described by the 
time series Xi and yi and is defined by the cross-recurrence 
matrix CR: 

CR,,,=e{e,-\\x,~y,\\), t,j^l,...,N, (1) 

where N is the length of the time series x and || • || is a 
norm, and O is the Heaviside function (Marwan et al. 120071) . 
The adjustable parameter Si is called the recurrence point 
density. 

In the special case of x — y, the main diagonal CR; i 
consists of ones ( CRi i = l,i = 1, ... , N). For this reason 
it is called the line of identity. Applying a transformation of 
the timescale (e.g., a phase shift) to the second time series 
y, the line of identity becomes bowed or shifted (Marwan 
et al. 120021 Marwan & Kurths ,2005) . It is caUed the line 
of synchronisation (LOS), because one can use it to rescale 
the timescale of the second time series and determine the 
closest match (synchronisation) of the two time series. 

The extraction of the LOS with the CRP is performed by 
the algorithm proposed in Marwan et al. ( 120021) . available 
in the CRP Toolbox for Madab® . When plotting the matrix 
with i and j being the axes of a diagram (see Sect. [3]), the 
LOS can be understood as a construction from the lower left 
to the upper right corner of the CRP-pattern, using the distri- 
bution of recurrence points within a window of size dx x dy 
with its origin at the identified recurrence point (i, j). 



3 Results 

From the analysis of the Greenwich sunspot area data it was 
shown that the consecutive sign changes in the phase dif- 
ference between the two hemispheres occurred at the maxi- 
mum of the 16th cycle and then at the minimum between the 
cycles 19 and 20 (Zolotova et al. 120091 ). At the same times 
the latitudinal distribution of the sunspots (expressed by the 
magnetic equator, defined by the sunspot latitudes) changed 
its sign (Pulkkinen et al. 119991 ). It is notable that the phase 
difference between the northern and southern hemispheres 
is in anticorrelation with the latitude of the magnetic equa- 
tor 

The raw daily sunspot counts reconstructed from the ob- 
servations by Carrington, Sporer, Staudacher, Hamilton, and 
Gimingham are not as uniform as the Greenwich data. Stau- 
dacher made his drawings at rather irregular time intervals. 
Only sometimes, he carried out more that 20-30 observa- 
tions during a year Carrington and Sporer usually did not 
admit gaps longer than 10 days, but they did not report day- 
by-day observations either. 

Averaging a time series with observational gaps is 
a commonly encountered problem when historical solar 
data are used (Vaquero 120071 Usoskin 2003). Since raw 
daily sunspot counts vary strongly, but we are interested 
in quite smooth curves, we used simple sliding averages 
over 500 days (for the Carrington and Sporer data) and 
over 1000 days (for the Staudacher, Hamilton, and Giming- 
ham data). Smoothing with such wide windows influences 
the cycle minimum and maximum levels, and causes time 
shifts. Thus, a strict numerical estimation of the amplitude 
of the northern and southern hemispheres or absolute phases 
may be incorrect. However, we are interested in phase dif- 
ferences and the hemispheres' behaviour during their as- 
cending and descending phases of the solar cycle evolution 
(Zolotova et al.2Q09). 

Smoothed sunspot counts reconstructed from the obser- 
vations by Carrington and Sporer for the north (blue), south 
(red), and the entire solar disk (black) are shown in Fig.[T] 
The smoothing procedure shortens the time-series by half 
of the window length at the beginning and the end of the 
raw data. Therefore, the smoothed Carrington and Sporer 
sunspot counts do not overlap (Fig.[Tl), although Carrington 
finished observing by the end of March 1861 and Sporer 
started in January 1861. During these about three months 
Carrington registered systematically more sunspots in com- 
parison to Sporer. We suggest that the discrepancies can 
be explained by a difference in detecting small sunspots or 
pores on the solar disc. 

The Sporer and Greenwich data overlap during a period 
of about 19.5 years. Both data sets were averaged with the 
same 500-day window. The Greenwich data are given by 
dashed lines in Fig.[T] During the maximum of the 12th cy- 
cle, Sporer registered on average two more sunspots in each 
hemisphere than the Greenwich observers. The next cycle 
does not show this difference. Nevertheless, both Sporer and 
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Fig. 1 Smoothed daily sunspot counts according to the observations by Carrington and Sporer for the north (blue), south 
(red), and the entire solar disk (black). The Greenwich data are marked by dashed lines in light colours (light-blue, pink 
and grey correspondingly). The cycle numbers refer to the Zurich numbering. 



Greenwich datasets show that from the beginning of cycle 
12, the northern hemisphere leads in time with respect to the 
southern. During the ascending phase of cycle 11, there is 
no hemisphere clearly leading, but the southern hemisphere 
shghtly predominates in leading continuously. The ascend- 
ing phase of cycle 10 also shows a leading southern hemi- 
sphere. Thus, we can claim the change of southern phase- 
leading into northern between the maximum of cycle 1 1 and 
the following minimum. We also confirm an earlier finding 
(Pulkkinen et al. 1 19991 Zolotova et al. |20091 l that the mag- 
netic equator lies in the hemisphere which is not leading in 
phase. 

In Fig. raw sunspot counts reconstructed from the 
observations by Staudacher, Hamilton, and Gimingham are 
shown for the northern (blue) and southern (red) hemi- 
spheres. Day-to-day observations (without gaps) are joined 
by lines. Some periods (solar maxima) are covered by fre- 
quent observations, but there are several gaps longer than 
half a year during solar minima. The average was not cal- 
culated if within the 1000-day window more than 990 days 
were without observations. Such periods correspond to the 
minima before the cycles 1 and 3, and the descending phase 
of the 4th cycle (Fig.|2}3, the entire solar disk sunspot counts 
ai-e shown in black colour) (Arlt l2008ll2009al l. 

Similarly to Zolotova & Ponyavin ('2006') and Zolotova 
et al. ([2009), to study the time-varying phase difference of 
the northern and southern hemispheres, we used the cross- 



recurrence plot (CRP) technique (Marwan et al. l2007l l. The 
algorithm is available in the CRP Toolbox for Matlab®. 

Since the original Staudacher, Hamilton, and Giming- 
ham records are too long for the CRP analysis, the data 
were resamplecQ. We extracted each 30th value from the 
smoothed daily sunspot counts. In a simplified attempt to 
construct the CRP which is sufficient for our purposes, 
these resampled numbers go into the time series Xi and 
Ui. Short resampled time series are practically identical to 
non-resampled, smoothed data shown in Fig.|2}3. Figure |2]; 
shows the cross-recurrence plot (gray pattern) of the resam- 
pled, smoothed sunspot counts with a diagonal line for ori- 
entation. The CRP was calculated using a fixed recurrence 
point density RR, e = 0.1. For the construction of LOS, the 
parameters dx = 30 and dy = 30 were used. 

LOS detects differences in timescales between two data 
series or dynamical systems (black thick line in Fig. |2j;). 
LOS variations relative to the zero level (or line of iden- 
tity) are shown in Fig. |2}l. The change in sign of LOS in- 
dicates the change of phase-leading from one hemisphere 
to the other (LOS > - north precedes south, LOS < - 

' Matlab (and the CRP toolbox as its application) operates slowly with 
loops for big matrices. The cross-recuirence plot is an N X N matrix 
(where A'^ is the time-series length). When N is more then some value 
(depends on the processor and RAM of the computer) Matlab toolbox 
stops operation or produces errors. A personal computer is able to process 
data with N equal to about 3500-5000. For the non-resampled Staudacher, 
Hamilton and Gimingham data. A' = 16826. Fortunately, the plots of non- 
resampled and resampled data ai'e highly similar. 
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Fig. 2 fl) Raw daily sunspot counts according to the ob- 
servations by Staudacher, Hamilton, and Gimingham for 
the northern (blue) and southern (red) hemispheres, b) 
Smoothed sunspot counts based on the same data set for the 
north (blue), south (red), and the entire solar disk (black). 
The cycle numbers refer to the Zurich numbering, c) Cross- 
recurrence plot of the resampled, smoothed sunspot counts. 
The Une of synchronization (LOS) is presented as black 
thick line, d) Extracted LOS. e) Latitude difference of 
sunspots distribution (the magnetic equator) versus time 
variations. Error bars indicated by gray colour. The vertical 
^sked/ilmevdcnfltgSntihachangevBfiilaflding near 1783 ±1.5 
yr. 



south precedes north). During the period from near 1750 to 
1783 ±1.5 yr (the uncertainty comes from half of the av- 
eraging window), the southern hemisphere predominantly 
leads. Afterwards, the entire 4th cycle is phase-led by 
the northern hemisphere. To account for the effect of fi- 
nite statistics, we used a crude approximation of errors 
^ l/C^^^, where C is the sunspot counts on the entire vis- 
ible solar disk during a day. 

Next, the magnetic equator was calculated as the sum of 
the mean latitudes (A) of sunspots in the northern and south- 
ern hemispheres {X{N))n + (A(5))„ (southern latitudes be- 
ing negative), where n is the time window, and is shown in 
Fig. 12^. The window size is the same - 1000 days, the re- 
quirement on free days from observations is also the same - 
no more than 990 out of 1000 days. The uncertainty in the 
sunspot latitude definition does not exceed 15° (Arlt l2009al l. 
It is seen that from about 1750 to 1783 the magnetic equa- 
tor was predominantly located in the northern hemisphere, 
while after that, in the beginning of cycle 4, it rapidly moved 
to the southern hemisphere. This again is in precise anti- 
correlation with the hemisphere leading in phase (vertical 
dashed line). 

Recently it was suggested that the unusual length of the 
4th solar cycle can be explained by an outstanding phase 
asynchrony between the northern and southern hemispheric 
activities reaching a delay of up to 4.5 years during the 
course of its ascending phase (Zolotova & Ponvavin l2007l l. 
However, from Fig. |2}l it is seen that the delay, during the 
first half of the cycle 4, between north and south does not 
exceed two years. The hypothesis about the lost cycle 4 
just before the Dalton minimum is actually discussed by 
Usoskin et al. (l2009l. 

During the Maunder minimum (1645-1715) the but- 
terfly symmetry was broken, and a significant ampli- 
tude difference between the hemispheres became evident. 
Sunspots were commonly observed in the southern hemi- 
sphere (Ribes & Nesme-Ribes 119931 Sokoloff & Nesme- 
Ribes'1994). But what about the phase difference of hemi- 
spheres during that epoch? 

Figure [3] is a copy of the butterfly diagram from Ribes 
& Nesme-Ribes (I1993I I. Vertical blue lines mark activity 
onsets and ends in the northern hemisphere, red lines - in 
the southern hemisphere. From 1683 to 1718 Philippe de 
la Hire conducted solar observations at the Observatoire de 
Paris. He carried out at least 10 observations (usually more) 
in each month. 

The north-south phase differences before the 1680s are 
difficult to fix, because sunspots lie very near the equator. 
In addition to that, the observations, before the 1680s, were 
irregular and rare (Ribes & Nesme-Ribes |19931 l. 

Near 1701 the activity in the southern hemisphere starts 
about two years earlier than in the opposite one. The next 
solar cycle starts near 1715 and the north leads the south 
by about one year. Thus, the sign change in the hemispheric 
phase difference probably occurred between 1 706 and 1715. 
However, up to 1714 the northern activity was very weak. 
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Fig. 3 Butterfly diagram of sunspots reconstructed for the period 1670-1719 according to Ribes & Nesme-Ribes ( |1993l l. 
Vertical coloured lines mark activity onsets and ends in hemispheres: north - blue, south - red. 



Hence, our hypothesis of the change of phase-leading from 
south to north should be considered a tentative suggestion. 
At the same time, up to 1714 the latitudinal distribution of 
sunspots was evidently shifted to the southern hemisphere. 

After the Maunder Minimum during the period from 
1715 to 1719 when the solar cycle recovered, the mag- 
netic equator cannot be determined easily simply by the eye. 
However, the northern activity precedes the southern one by 
about one year. 

4 Discussion 

We investigate variations of the hemispheric phase differ- 
ence in the solar cycle in the period of 1670-2009 by a com- 
bination of positional sunspot data from various sources. 
The results are shown in Fig. |4] together with the yearly 
group sunspot numbers by Hoyt & Schatten (1 19981 1 for 
1610-1995. 

We can summarize that during 1715-1719 the northern- 
hemisphere cycle precedes the southern one (Fig. H] corre- 
sponding to the time interval coloured by light-blue). Since 
this conclusion was made only from a visual analysis of the 
butterfly diagram during the Maunder minimum, it is not 
a robust result. The northern-hemisphere leading probably 
lasts 30^5 years until a sign change in the phase differ- 
ence occurs near 1745-1750. Then, according to the Stau- 
daucher's drawings, the southern hemisphere precedes for 
35-40 years and hands over to the northern hemisphere 
near 1783. The next recurrent change in phase-leading is 
supposed to lie somewhere between the maxima of the 7th 
and 8th cycles (grey bar near 1833). This change in phase- 
leading is not documented by data yet. Such a long interval 
for a sign change (roughly 10 years) is chosen because the 
secular variations of phase asymmetry as well as of the mag- 
netic equator are not symmetric sinusoidal functions. An in- 
dividual phase-leading period by either hemisphere may last 



somewhat longer or somewhat shorter than 4 cycles. Thus, 
the northern hemisphere preceded the southern one from the 
beginning of the cycle 4, probably for 45-55 years. Accord- 
ingly, the southern hemisphere then preceded the northern 
one for 38^8 years. 

According to the Carrington and Sporer data, the south- 
ern hemisphere leads in phase during the 10th cycle and 
during the ascending phase of the 1 1th cycle, but by the on- 
set of the 12th cycle the northern hemisphere starts to pre- 
cede the southern one. This sign-change in phase difference 
happened gradually, during the descending part of the 1 1th 
cycle. Note that, in contrast, another north-south change in 
phase-leading occurring in the minimum between cycles 19 
and 20 was abrupt (Zolotova et al. ,2009j . 

Further long-term variations in the hemispheric phase 
differences are described by Zolotova et al. ( 120091 1: from 
the beginning of the 12th to the maximum of the 16th cy- 
cle, the northern hemisphere dominates in leading for about 
50-55 years period, after which, until the minimum before 
the cycle 20, the southern hemisphere leads for about 35-40 
years, and then, till the present, in spite of depressed sunspot 
activity during the new solar cycle 24, the northern hemi- 
sphere leads in time. To keep pace with the proposed secu- 
lar variation in hemispheric leading, we expect the next sign 
change in phase difference as well as the northward shift of 
the magnetic equator will happen no later than in the end of 
the current solar cycle (2018-2020). For convenience, the 
results are also shown in Table 1 . The hemispheres are not 
equal in their duration of phase-leading. The duration of the 
northern hemisphere preceding the southern one is longer 
than the opposite case. The cause of such a skewness is un- 
clear. 

From 1750 to 1796 and from 1853 to the present, phase- 
leading is found to be antic orrelated with the latitudinal 
distribution of the sunspots. That is, the magnetic equator 
will be located in the hemisphere which lags behind the 
other hemisphere in its activity. However, such relationship 



www.an-journal.org 



© XXX WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 



6 



N.V. Zolotova et al.: Secular variation of liemispheric pliase differences in the solar cycle 



35-40 
years 



45-55 
years 



38^8 
years . 



35^0 45-50 
years I years? 




1600 1650 1700 1750 



1800 
Year 



1850 1900 1950 2000 



Fig. 4 Yearly group sunspot numbers versus time. Periods of the northern hemisphere preceding the southern one are 
marked by blue colour, periods when the southern hemisphere precedes the northern one are marked by red colour. Vertical 
lines denote the sign changes of the phase difference. Grey bars indicate probable sign changes. The cycle numbers refer 
to the Zurich numbering. 



Table 1 Preceding hemispheres over the last 300 years. 



Preceding Onset End Duration 

Hemisphere [yr] 



Northern 


1706- 


-1715 


1745- 


-1750 


30-45 


Southern 


1745- 


4750 


1783 


35-40 


Northern 


1783 


1829- 


-1837 


45-55 


Southern 


1829- 


4837 


1873- 


-1877 


38^8 


Northern 


1873- 


-1877 


1928 


50-55 


Southern 


1928 


1968 


35-40 


Northern 


1968 









is not continuous, and was broken during the Maunder min- 
imum. Probably it can be explained by qualitatively differ- 
ent dynamo behaviour during Grand minima compared with 
the normal dynamo-states (Brandenburg & Spiegel 120081 
Brooke et al. lT998l Kiiker et al. 119991 Ossendrijver|2003ll. 

Note that, interpreting the solar cycle as a pair of activity 
waves travelling towards the equator and pole respectively. 
Pelt et al. (I2000I I found the Gleissberg cycle modulation not 
only in the variations of the magnetic equator, but in the 
period length of the Schwabe cycle as well. 



5 Conclusion 

We used historical documents by Staudacher, Hamilton, 
Gimingham, Carrington, Sporer, and the Greenwich ob- 
servers, as well as the sunspot activity reconstruction by 
Ribes & Nesme-Ribes ( |19931 l. We showed that one of the 
two hemispheres is always shghtly advanced in the progress 
of its cyclic activity compared to the other hemisphere dur- 
ing the last 300 years. After several cycles, the hemispheres 
change their roles. The long-term persistence of this phase- 
leading does not demonstrate a clear periodicity, but varies 
in duration around 8 solar cycles. The periodicity proba- 
bly coiTesponds to the Gleissberg cycle. Phase-leading is 
found to be antic orrelated with the latitudinal distribution 
of the sunspots. This means that the sunspots in a preceding 
hemisphere show a butterfly wing emerging on a slightly 
lower latitudes as compared to the butterfly wing of the op- 
posite, delayed hemisphere, where the sunspots emerged at 
slightly higher latitudes. Since the north-south asymmetry 
exhibits a long-term persistence, it should not be regarded 
as a stochastic phenomenon. 

Acknowledgements. For the calculation of CRP and LOS, 
we have used the CRP Toolbox, developed by Norbert 
Marwan, for Matlab®, available at tocsy.pik-potsdam.de. 
The Royal Greenwich Observatory USAF/NOAA daily 
data of sunspot area and their latitudes separately for 
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the northern and southern hemispheres is available at 
http://solarscience.msfc.nasa.gov/greenwch.shtml The group 
sunspot numbers (GSN) data provided by Hoyt & Schat- 
ten ( |1998l l and tabulated by the National Geophysical Data 
Center is available at http ://www.ngdc.noaa.gov/stp/SOLAR/ 
ftpsunspotnumber.html. 
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